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ABSTRACT 

Quantitative structural analysis of the galaxies present in the Hawaiian Deep Fields 
SSA13 and SSA22 is reported. The structural parameters of the galaxies have been 
obtained automatically by fitting a two-component model (Sersic r 1 /™ bulge and ex- 
ponential disc) to the surface brightness of the galaxies. The galaxies were classified on 
the basis of the bulge-to-total luminosity ratio (B /T) . The magnitude selection criteria 
and the reliability of our method have been checked by using Monte Carlo simulations. 
A complete sample of objects up to redshift 0.8 has been achieved. Spheroidal objects 
(E/SO) represent » 33% and spirals w 41% of the total number of galaxies, while 
mergers and unclassified objects represent 26%. We have computed the comoving 
space density of the different kinds of objects. In an Einstein-de Sitter universe a 
decrease in the comoving density of E/S0 galaxies is observed as redshift increases (a 
w 30% less at z=0.8), while for spiral galaxies a relatively quiet evolution is reported. 
The framework of hierarchical clustering evolution models of galaxies seems to be the 
most appropriate to explain our results. 

Key words: galaxies: distances and redshift — galaxies: evolution — galaxies: 
photometry — galaxies: fundamental parameters 



1 INTRODUCTION 

Achieving a good galactic evolutionary model is one of the 
challenges of present astronomy. The high quality of the 
Hubble Space Telescope (HST) data allows astronomers to 
study the evolution of galaxy morphology over a signifi- 
cant fraction of the age of the Universe, restricting the two 
main present theoretical frameworks of galaxy evolution: the 
monolithic collapse and the hierarchical clustering models. 

The simplest models of galaxy evolution predict that 
massive elliptical galaxies formed at high redshift in a rapid 
collapse with a single burst of star formation (Eggen et al 
1962; Larson 1975). Against this scenario, the hierarchical 
clustering models predict that the most massive objects form 
at late times via the merging of smaller subunits (White & 
Rees 1978; Kauffmann et al 1993). Each model has very 
different observational implications (e.g. Brinchmann et al 
1998; Schade et al 1999; Fried et al 2001). Observational ev- 
idence has been found for both scenarios (see Schade et al 
1999 and references therein), so that the dominant mecha- 
nism of galaxy evolution remains an open question. 

Many attempts have been made to classify galaxies on 
HST deep images. Two families of methods are currently 
used: visual and automated classifications. Among visual 
classifications we mention analysis done by van der Bergh 



et al (1996, 2000) in the range 21<I 8 i 4 <25 at the Hubble 
Deep Field (HDF). They found that up to 30% of the galax- 
ies were ellipticals, the remainder divided into 31% spirals 
and 39% unclassified. Possible differences in the morpholo- 
gies of galaxies at high redshifts point to different the envi- 
ronmental conditions of these galaxies relative to the local 
ones. In particular, the merger rate could be very different. 
Le Fevre et al (1999) have found that the rate of mergers 
and interaction grows strongly with the redshift. Quantita- 
tive classification systems based on the study of the central 
concentration and asymmetry of the galaxian light (Abra- 
ham et al 1996) also obtained a high fraction of irregular 
and peculiar galaxies at high redshifts, finding only a 20% 
of elliptical fraction. 

Most sophisticated classification systems based on the 
decomposition of the surface brightness profiles of galaxies 
into their structural components (bulge and disc mainly) 
have been applied during the last few years. This technique 
is used extensively for local galaxies (see Prieto et al 2001 
and references therein) but the lower resolution at high red- 
shift makes its application there more difficult . This quan- 
titative classification method has the advantage that it gives 
information about each component of galaxies. This means 
that we can follow the evolution of different components 
(bulge and disc in spirals) as a function of redshift. Usu- 
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ally, it is assumed that the same type of profiles which fit 
the light distribution of local galaxies also describe the light 
distribution of galaxies at higher redshift. Typically, Sersic 
j- 1 /™ profiles are fitted to the surface brightness profiles of 
bulges and elliptical galaxies, and exponential profiles to the 
discs of Spiral galaxies (Marleau & Simard 1998, Schade et 
al 1996, 1999). 

Using this decomposition technique on the HDF, Mar- 
leau & Simard (1998) found a substantially different result 
from those obtained by visual classifications. They found 
that only 8% (versus 30% for visual classifications) of the 
galaxy population down to Isu{AB) — 26 are spheroidal 
systems. Although quantitative methods have clear advan- 
tages over visual methods, they are not free from significant 
bias which affect the reliability of the physical properties ob- 
tained. In order to understand the big discrepancy pointing 
out in the previous analysis, it is crucial to remove the bias 
which are present in quantitative analysis methods. 

In this paper, we examine the structural properties 
of the galaxies in two Hawaiian Deep Fields (SSA13 and 
SSA22) imaged by HST. Each of these fields is composed 
by 3 HST/WFPC2 fields. All the galaxies studied in these 
fields have spectroscopic redshifts, avoiding a strong source 
of uncertainty at the distance determination. Previous clas- 
sification schemes of high redshift galaxies from HST images 
are compared with our results. In particular, we focus our 
attention on evaluating the number of spheroidal systems in 
field galaxies and on constraining the two main theoretical 
frameworks of galaxy evolution. 

The structure of this paper is as follow: section 2 de- 
scribes the characteristics of the observed fields. The struc- 
tural decomposition method is presented in section 3. In 
section 4 we discuss the completeness of the sample and we 
summarize our conclusions and discuss their implications for 
galaxy evolution in sections 5 and 6. 
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2 THE SAMPLE OF GALAXIES 

The sample consists on all objects with K<20, I<22.5 (Kron- 
Cousins) and B<24.5 in two areas surrounding the Hawaii 
deep survey fields SSA 13 and SSA 22 (Cowie et al 1994; 
Songaila et al 1994). Hereafter, the I magnitude will be given 
in the same system as in Cowie et al (1996). Nearly all ob- 
jects included in those fields have measured spectra with the 
LRIS spectrograph on Keck (Cowie et al 1996) The fields 
were imaged during 2000 s with the WFPC2 at HST in the 
1814 bandpass. The total analyzed sky area was 28 arcmin 2 . 
The analyzed objects lie in the redshift interval [0.1, 1.3], 
mainly concentrated around z = 0.5 (see Cowie et al 1996). 

We used the SExtractor galaxy photometry package 
(Bertin & Arnouts 1996, version 2.1.4) for the extraction 
of the objects from the public released HST images. This 
package is optimized to detect and measure sources from 
astronomical images. The detection was run using the same 
parameters as in Marleau & Simard (1998). In particular, we 
used a detection threshold of 1.5 a, where a is the standard 
deviation of the sky background of the images. Another im- 
portant parameter is the deblending parameter. SExtractor 

* See the discussion about the different magnitudes and trans- 
formations in Cowie et al (1995) 



mag(input) 

Figure 1. (Top) d(S/T)=B/T(measured)- J B/T(input) as a func- 
tion of the input magnitude. (Bottom) Mean d(B/T) vs. input 
magnitude with 1 a error bars. 



deblends objects using multiple flux thresholding. The SEx- 
tractor deblending parameter sets the minimum fraction of 
the total flux a branch must contain to be considered a sep- 
arate object. We have use the same value that Marleau & 
Simard (1998), which is 0.001. 

In order to obtain a bulge+disc decomposition of the 
objects, we fit ellipses to their isoluminosity contours down 
to 1.5 a using the task ELLIPSE from IRAF. The surface 
brightness and ellipticity profiles obtained are used to re- 
cover the structural parameters of the galaxies. 



3 THE GALAXY CLASSIFICATION 
PROCEDURE 

The classification technique is based on the decomposition 
of the surface brightness profiles of the galaxies in bulge and 
disc components. The fitting algorithm is discussed exten- 
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sively in Trujillo et al (2001b). Here we explain the main 
points ol the routine. 

The final surface brightness distributions resulting from 
the convolution between the PSF and our 2D (i.e. elliptical) 
model surface brightness distributions are dependent on the 
intrinsic ellipticity of the original source - as is the case with 
real data. A key problem remains, which is what value of the 
ellipticity is choosen to represent the ellipticity of the model. 
The ellipticity of the isophotes are reduced by seeing. This 
reduction depends on the radial distance of the isophote to 
the center of the model, the size of the seeing, and the values 
of the model parameters. Consequently, to evaluate the in- 
trinsic ellipticity of a model it is often insufficient to simply 
measure the ellipticity at one given radial distance (e.g. 2 
effective radii). To illustrate this, the observed ellipticity at 
2 r e on galaxies which have an effective radius of similar size 
to the FWHM (these galaxies are common at high redshift) 
is 30% less than the true ellipticity for galaxies with an ex- 
ponential profile (n=l), and 45% less for galaxies with a de 
Vaucouleurs profile (n=4). The use of models with underes- 
timated ellipticity affects the evaluation of the other model 
parameters, biasing the results. One result of this bias is the 
estimation of smaller values of index n. This bias increases 
as the value of n increases (Trujillo et al. 2001a,c). 

Consequently, the determination of the intrinsic ellip- 
ticity of the source and the fitting process to determine 
the structural parameters should be done in tandem (i.e. 
using an iterative and self-consistent routine) and not as 
two separate tasks. To do this we simultaneously fit both 
the observed surface brightness and ellipticity profiles using 
convolved profiles for each (see how the algorithm works in 
Figure 6 from Trujillo et al. 2001b). 

Our 2D fitted galaxy model has two components: a 
bulge and a disc. The 2D bulge component is a pure Sersic 
(1968) profile of the form[} 

-M(€A„) 1/ "-i] 



7(C) = /elO" 



(1) 



where I e is the effective intensity, r e is the bulge effec- 
tive radius and b n = 0.868n - 0.142 (Caon et al 1993). The 
disc component is an exponential profile given by: 

1(0 = I e^ h (2) 

where I is the central intensity and h is the exponential 
disc scale-length. The set of free parameters is completed 
with the ellipticities of the bulge e;, and the disc ed- The 
bulge and disc profiles were convolved with the instrumen- 
tal PSF of the HST obtained from stellar profiles located on 
the images. Special attention was paid to this convolution. 
The real PSFs were fitted by Moffat functions and the con- 
volutions were developed analytically on real space. Also, to 
avoid the problem of the undersampling of the PSF we av- 
erage different stellar profiles obtaining a composed median 
PSF. To this median profile we fit our analytical PSF. We 
have estimated a ~5% uncertainty in the estimation of the 
FWHM due to changes from one WFPC2 position to an- 
other. This uncertainty implies an error on the parameters 
estimation less than 10%. 

A Levenberg-Marquardt non-linear fitting algorithm 

t The surface brightness distribution are explicitily written on 
elliptical coordinates (£, 8) (Trujillo et al. 2001a). 
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Figure 2. (Top) The difference, d(B/T), between measured and 
input B/T versus B/T(input) for two different magnitude inter- 
vals: I < 23 (solid circles) and 1^23 (triangles). (Bottom) Mean 
d(B/T) versus B/T(input) with 1 a error bars. 



(Press et al 1992) was used to determine the free parameters 
set which minimizes \ 2 - Extensive Monte-Carlo simulations 
were done in order to check the reliability of the recovered 
parameters (see Section 4). The surface brightness profiles 
and ellipticity profiles of each galaxy were fitted at the same 
time. Each galaxy was fitted by a single Sersic profile and a 
Sersic + exponential profile. 

Following previous studies (e.g. Marleau & Simard 
1998), galaxy classification was based on the bulge to to- 
tal luminosity ratio, B/T. We consider as "ellipticals" those 
objects with B/T > 0.6, in which case a better fit can be 
obtained with only one component. The parameters of these 
objects were taken from the pure Sersic fitting. Galaxies with 
B/T between 0.5-0.6 were classified as SO. Finally, objects 
with B/T < 0.5 were classified as "spirals". We consider 
as "spheroidal" galaxies those with B/T >0.5 as Marleau 
& Simard (1998). The discrimination between the different 
types of galaxies was made following the values of B/T given 
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by Simien & de Vaucouleurs (1986). Quantitative selected 
ellipticals can be contaminated by galaxies such as com- 
pact narrow emission-line objects. These galaxies exhibit 
high bulge fractions even though they are not real ellipticals. 
These objects can be ~15% of the elliptical sample (Im et al. 
2001) and their presence must be taken into account at es- 
timating the uncertainty on the ellipticals comoving density 
parameters^]. 

Galaxies selected using only the B/T >0.5 criteria may 
not all be E/SOs, but could include later galaxy types. To 
quantify this bias we use the analysis performed by Im et al. 
(2001) for a local galaxy sample (Frei et al. 1999). Galaxies 
with T <0 (i.e. E/SOs) represent 76% of the local sample 
selected using B/T >0.5. So, a contamination of ~ 25% can 
be expected in the objects that we are labeling as E/SOs at 
high redshift. However, the contamination for objects with 
T <0 in the objects named "spirals" (i.e. B/T <0.5) is just 
8%. Some methods have been identified to remove the bias in 
the E/SOs selected sample with the use of red colors selected 
galaxies or the use of low asymmetry objects. However, the 
first option clearly biases the sample to objects that have 
a quiet evolution (and what we want is precisely study this 
hypothesis) and the second has been shown to be inappro- 
priate in objects at high redshift (i.e. low S/N as our objects 
have) by Conselice, Bershady & Jangren (2000). Despite the 
known morphological type biases, due to the above reasons 
we have chosen to maintain the B/T selection criteria as the 
sole morphological selection criteria. 

For the ellipticals, we have also imposed a restriction 
based on absolute magnitude. By doing this, we have clas- 
sified a galaxy as "dwarf" when Mb > —17.0. The absolute 
magnitudes were obtained after applying the K-correction 
prescription of Poggianti (1997) assuming (hereafter) a cos- 
mology with Ho=75 km s _1 Mpc -1 , qo=0.5, Q m = 1.0 and 
O\=0. 

Once the automated classification is done, a visual in- 
spection was also made for each object. Some objects are not 
fitted well by either a pure Sersic profile or a bulge+disc pro- 
file. They were classified as "irregular" galaxies. Those with 
evidence of mergers (close companions and irregular shapes) 
were catalogued as "mergers". We also had 4 objects whose 
best fit is achieved by a pure Sersic profile with n w 0.5. It 
is important to note that the luminosity density of a Sersic 
profile with n< 0.5 has a depression in its central part repre- 
senting an unlikely physical situation (Trujillo et al 2001a). 
Marleau & Simard (1998) also obtained some objects of this 
class on the HDF images. The visual morphological shape 
of these objects is peculiar appearing elongated. Marleau & 
Simard (1998) claimed that this kind of objects could be 
remnants of mergers or close tidal disruptions. We have in- 
cluded them into the merger category. 
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T As a matter of caution we also must regard that these "inter- 
lopers" are basically placed at high redshift (z>0.8) or they are 
faint galaxies Mb <-18 (see Im et al. 2001 their Figure 17). For 
that reason, most of these galaxies are expected to be out of our 
studied sample. 



Figure 3. The complete magnitude as a function of the redshift 
for objects with apparent magnitudes: I = 22 (full line) and / = 
23 (dashed line). Three different kind of objects are represented: 
ellipticals (top), Sa (middle) and Sc (bottom). See text for details. 



4 THE COMPLETENESS OF THE SAMPLE 

Since selection effects can mimic evolutionary changes in 
high redshift objects, it is necessary to achieve a complete 
unbiased sample of objects. The determination of the com- 
pleteness of the sample is done in two steps. First, we de- 
termine the faintest apparent magnitude down to which the 
recovered parameters are reliable. In particular, we will fo- 
cus on the B/T ratio because it is the parameter used for 
the classification of the galaxies. We evaluated this limiting 
magnitude by Monte-Carlo simulations of artificial galaxies 
with similar magnitudes and structural parameters as the 
real objects. Once this magnitude is obtained, the second 
step for the completeness of the sample consists in deter- 
mining how bright (i.e. the absolute magnitude) a galaxy 
has to be in order to be observed in our whole redshift inter- 
val. The limiting absolute magnitude was obtained by using 
typical spectra from every type of objects, which allows us 
to verify that we are studying the same kind of objects in 
all redshift intervals. Unfortunately, most previous studies 
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Figure 4. The Mi — z diagram for the galaxies detected at the 
SSA13 and SSA22 fields. The absolute magnitudes have been 
computed for a H = 75kms~ 1 Mpc~ 1 , f2 m = 1 and q a = 0.5 
cosmology. 



of the structural properties of high redshift samples do not 
determine their limiting absolute magnitudes. This type of 
samples are obtained with only an apparent limiting mag- 
nitude, which biases the sample to the brightest objects at 
high redshifts. For this reason, it is necessary to use a sample 
cut by absolute magnitude. 

4.1 Monte-Carlo simulations 

We performed Monte-Carlo simulations to test the reliabil- 
ity of our method. First, we tested the ability to recover pa- 
rameters from bulge-only (i.e. purely elliptical) structures, 
and second we explored the possibility of carrying out ac- 
curate bulge+disc decompositions. In both cases we created 
150 artificial galaxies with structural parameters randomly 
distributed in the following ranges: 

• bulge-only structures: 19<I<23, 0.05" <r e <0.6", 0.5< 
n < 4, and < e < 0.6 (the lower limit on n is due to the 
physical restrictions pointed out in Trujillo et al 2001a). 

• bulge+disc structures: 18.5<I<23.5, 0.05"<r e <0.6", 
0.5< n < 4, and < e b < 0.4, 0.2" < h < 1.5", 0< B/T < 1, 
and < e d < 0.6. 

The artificial galaxies were created by using the IRAF 
task MKOBJECT. We support as an input to this task 
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Figure 5. The cumulative number distribution function of E/SOs 
(top) and spirals (bottom) as a function of the redshift. The best 
fits derived from linear and quadratic comoving densities are over- 
plotted. It is also showed the cumulative distributions derived 
from Fried et al. (2001). See text for more details. 



the surface brightness distribution coming from our detailed 
convolution between the PSF and the original model. To 
simulate the real conditions of our observations, we added a 
background sky image (free of sources) taken from a piece 
of the real image; the dispersion in the sky determination 
was 0.1 %. The PSF FWHM in the simulation was set at 
0.2" and assumed known exactly. The pixel scale of the sim- 
ulation was 0.1", as is the real WFPC2 pixel size. The same 
procedure was used to process both the simulated and the 
actual data. 

Figure 1 shows d(B/T)=B/T(measured)-B/r(input) 
as a function of the input magnitude. For galaxies brighter 
than I = 23 magnitude d(B /T) is less than 0.1. This is 
a very accurate determination of this parameter. Figure 2 
shows d(B/T) as a function of the B /T(input). Objects with 
I > 23 (triangles) have bigger dispersion of d(B/T) than ob- 
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Table 1. Galaxies from SSA13 with / < 23.0. 



Table 2. Galaxies from SSA22 with / < 23.0. 
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jects with / < 23 (points). From our simulations it follows 
that / = 23 is the limiting magnitude for reliable recover 
of the B/T parameter. The limiting magnitude for the rest 
of structural parameters will be studied in a forth-coming 
paper. To our limiting apparent magnitude, the sample of 
galaxies is reduced to 120 galaxies. According to their B/T 
ratio, absolute magnitude and visual inspection (see Section 
3), they were classified as: ellipticals (26), dwarfs (6), SO 
(9), irregulars (20), mergers (17) and spirals (42). This left 
us with ~ 34% spheroidal galaxies (E+Dwarfs+S0), ~ 35% 
spirals and ~ 31% of unclassified objects. Tables 1 and 2 
show the B/T ratios (column 4) and Mi (column 3) for the 
120 galaxies with I < 23.0 for the SSA13 and SSA22, re- 
spectively. These tables also show the identification number 
(column 1) and the redshift of the objects (column 2) given 
by Cowie et al (1996). Galaxies classified as ellipticals and 
dwarfs have B/T = 1.0, those classified as irregular of merg- 
ers have been marked with B/T = —1.0. 

4.2 Completeness as function of redshift 

In order to be sure that we are studying the same kind of 
objects at different redshifts we must determine the absolute 
limiting magnitude of our sample. On doing this we avoid 
biasing our sample to brighter objects at high redshift. Some 
claims of galactic evolution have been a consequence of this 
bias. As an example, Simard et al (1999) analyzed the prob- 
lem of the completeness of the sample. If selection effects 
were ignored in their galaxies, then the mean disc surface 
brightness increases by ~ 1.3 magnitudes from z = 0.1 to 
z = 0.9. Most of this evolution is plausibly due to compar- 
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-19.58 


-1.00 


89 


1.151 


-22.46 


0.41 


28 


0.247 


-21.41 


0.04 


90 


0.412 


-19.52 


1.00 


30 


0.751 


-21.78 


0.34 


91 


0.513 


-20.13 


1.00 


32 


1.024 


-22.54 


-1.00 


92 


0.381 


-20.94 


0.50 


33 


0.707 


-22.07 


1.00 


93 


0.377 


-19.95 


1.00 


38 


1.208 


-24.05 


1.00 


96 


0.290 


-22.21 


1.00 


41 


0.672 


-21.88 


0.51 


100 


0.303 


-19.93 


1.00 


45 


0.132 


-18.15 


0.53 


102 


0.824 


-22.39 


0.20 


46 


0.912 


-22.04 


0.08 


103 


1.159 


-23.70 


1.00 


47 


0.173 


-17.38 


1.00 


108 


0.588 


-21.25 


-1.00 


19 


0.707 


-20.98 


0.37 


111 


0.302 


-18.47 


0.21 


50 


0.538 


-21.22 


0.09 


118 


0.816 


-21.36 


0.36 


51 


0.536 


-21.11 


-1.00 


123 


0.095 


-19.04 


0.11 


54 


0.418 


-20.67 


-1.00 


124 


0.671 


-20.81 


-1.00 


55 


0.815 


-21.58 


-1.00 


125 


0.873 


-22.00 


0.24 


56 


0.318 


-17.80 


0.23 


127 


0.695 


-21.39 


1.00 


59 


0.418 


-20.58 


0.50 


143 


1.102 


-22.31 


0.09 


60 


1.392 


-19.54 


-1.00 


147 


0.514 


-21.24 


1.00 


67 


0.588 


-19.45 


-1.00 


148 


0.876 


-22.41 


-1.00 


69 


0.692 


-21.67 


1.00 


150 


0.795 


-22.11 


1.00 


70 


0.348 


-20.62 


0.38 


152 


0.617 


-22.27 


0.17 


71 


0.132 


-16.72 


1.00 


154 


0.614 


-22.09 


1.00 


72 


0.787 


-21.07 


0.07 


155 


0.665 


-20.71 


1.00 


73 


0.822 


-23.08 


1.00 


161 


0.960 


-22.65 


1.00 


75 


0.724 


-20.13 


-1.00 


166 


0.378 


-18.76 


1.00 


77 


1.020 


-22.11 


0.36 


172 


0.378 


-19.93 


-1.00 


78 


0.823 


-21.49 


0.55 


204 


0.709 


-19.77 


-1.00 



ing low-luminosity galaxies in nearby redshift bins to high- 
luminosity galaxies in distant bins. If this effect is taken into 
account, no discernible evolution remains in the disc surface 
brightness of their disk dominated galaxies. In order to avoid 
this kind of problem it is necessary to make a selection of 
the objects on the basis of their absolute luminosity. 

Given our apparent limiting magnitude of / = 23, we 
have studied the completeness of our sample for three dif- 
ferent class of galaxies: Sa, Sc and E. Figure 3 shows, for a 
limiting magnitude of 1=23 , the absolute magnitude down 
to which a galaxy can be observed as a function of z. This 
figure was generated using spectral models of 15 Gyr old 
galaxies (Poggianti 1997). For our distribution of 120 ob- 
jects with I < 23, the parameters which maximize the num- 
ber of objects into a complete sample are z < 0.8 and 
Mi < -20.0 (M B < -18)|| This left us with a total of 
61 objects: 20 E/S0, 25 spirals, 16 irregulars and mergers, 
or equivalently:~ 33% E/S0 ~ 41% spirals and ~ 26% un- 
classified objects. Figure 4 shows the M — z diagram for 
our whole sample down to / < 23. The complete subsample 
studied (left-down corner) is enclosed by the horizontal and 
vertical lines. This kind of selection criteria is similar to that 
used by Simard et al (1999) and Fried et al (2001). 

§ We have repeated this calculation with the starburst galaxy 
NGC4449 without finding any substantial difference. 
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5 RESULTS AND DISCUSSION 
5.1 Galaxy classification 

There is no substantial differences between the fraction of 
galaxies which correspond to the different classes when the 
sample is restricted in apparent and absolute magnitude. We 
must note that these numbers are in good agreement with 
the percentage of E/SO given by visual classification systems 
(van den Bergh et al 1996, 2000) in the HDF, but are quite 
different from that given automated classification of Mar- 
leau & Simard (1998), 8%. The classification of Marleau & 
Simard (1998) takes into account all objects with 7si4 < 26 
from the HDF. They claim that the discrepancy with visual 
classifications is due to the difference in the classification 
of small round galaxies with half-light radii less than 0."31. 
Visually these galaxies are classified as elliptical galaxies, 
Marleau & Simard classify them as disc-dominated systems 
with bulge fractions less than 0.5. But, galaxies with an 
intrinsically big B/T tends to be systematically obtained 
with lower values of B/T in their automated routine (see 
Fig. 10 of Marleau & Simard 1998). The ellipticals in the 
HDF which have been (probably) mis-classified by these 
authors are those principally coming from the fainter sub- 
sample. Although the relation between the B/T output of 
their simulation and the input magnitude of the objects is 
not provided by these authors (and consequently our asser- 
tions must be taken with caution), it is certainly possible 
that B/T results stronger affected at increasing the input 
magnitude (i.e. at lower signal ratios), and for that reason, 
the high redshift population of elliptical galaxies remains 
biased. By using an automated procedure which avoids this 
problem we have been able to obtain a result similar to van 
den Bergh et al (1996, 2000). 

Interestingly, our sample and the HDF one are imaging 
a galaxy population centered around z ~ 0.5, the principal 
difference being in the exposure time. Because of the dif- 
ferent depth in the images, substantial differences would be 
expected for the fainter subpopulation (smaller and irreg- 
ular galaxies) between our sample and those based on the 
HDF. In fact, in the HDF apparent magnitude-limited sam- 
ple contains 39% of unclassified objects whereas we obtain 
- 30%. 



5.2 Galaxy evolution 

The two main models of galaxy evolution (monolithic col- 
lapse and hierarchical clustering) present a completely differ- 
ent scenario of galaxy evolution, so that the observational 
implications also are very different. One of these concerns 
the comoving density of the galaxies. In the redshift interval 
studied, the hierarchical model framework proposes the co- 
moving density of big galaxies (E/SOs and spirals) decrease 
with redshift, being constant in the monolithic model. We 
have computed the comoving density p(z) for these two dif- 
ferent types of galaxies in our complete subsample. 

We have assumed both a linear function for modeling 
the comoving density, 

p(z) = a + bz (3) 

with a=p(0) the comoving density at z~0 and 
b=(p(z max ) — p(0))/zmax with z max the maximum value 



which z reaches in our limiting subsample, and a power law 
of the form: 

p(z) = a(l + z) b (4) 

with a=p(0) the comoving density at z=0. 

To reduce the loss of information in our data we avoid 
binning them. The values of the parameters of the function 
p(z) are achieved by running a Kolmogorov-Smirnov test 
between the cumulative probability distribution function of 
finding a galaxy inside our imaging solid angle at a given 
z given p(z) and the cumulative distribution from the real 
data. The cumulative probability function for our model is 
computed as: 

P p(z')r 2 (z')(dr/dz)dz' 

p( z ) = Zm '" (5) 

J*™* p(z')r 2 (z')(dr/dz)dz' 

where z m i n is the closest galaxy redshift, z ma x=0.8 for our 
limiting subsample and r(z) is the comoving distance to an 
object placed at z. 

The Kolmogorov-Smirnov (KS) test gives the probabil- 
ity that two data sets come from the same distribution. The 
best comoving density is that which maximize this proba- 
bility. 

As a matter of caution, we must note that E/S0 galaxies 
are placed preferentially in high density environments, being 
more strongly clustered than other types of galaxies. In order 
to evaluate the effect of clustering in our comoving density, 
we have studied the contribution of the E/S0 galaxies of each 
field to our cumulative function. The number of galaxies 
on the regions of more accumulation (0.48-0.54 and 0.69- 
0.73) come from both fields with approximately the same 
contribution, rejecting a clustering explanation. 

To evaluate the errors on the parameters in the E/S0 
and in the spiral sample we have assumed that a least two 
galaxies in each sample are misclassified. This represents 
~10% of each sample. We construct all the subsamples that 
can be obtained by removing two elements on the original 
samples and then we recover the values of the parameters as- 
sociated to them. Using these values we estimate the median 
and the standard deviation. These are the numbers that we 
present as the parameter estimations and the errors associ- 
ated to these measurements. Fig. 5 shows the original whole 
sample (i.e. without removing any point) and overplotted 
is the cumulative function associated with the parameters 
measured as explained before. Bar errors in Fig. 5 were es- 
timated by measuring at each point the maximum distance 
between the cumulative function represented by using the 
whole sample and all the cumulative functions resulting from 
the previous subsamples. We have also overplotted the cu- 
mulative distribution obtained from the comoving densities 
fitted by Fried et al. (2001), who have a similar absolute 
magnitude cut for their sample (Mb < —18.5). 

The comoving density of the E/SOs which gives a max- 
imum probability in the KS test for a linear form is given 
by: p(z) = 0.0033(±0.0015) - 0.0015(±0.0010) x z. The KS 
probability of this density is 0.90. This comoving density 
is closer to that deduced by Fried et al (2001). Using their 
fit to our sample we obtain a KS p (KS pf) of 0.87. The 
number of E/SOs decreases with redshift. For the cosmology 
chosen, this decrease is ~ 45(±30)%. At using the power-law 
model, the KS probability is slightly better, 0.92. We have: 
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Table 3. Parameters of the comoving density. Ho=75kms — ^^Mpc - 1 



Model 


Type 


a 


b 


c 


KS p 


KS pf 






fi m — 1 




Q A = 






a+bxz 


E/SO 


0.0033(±0.0015) 


-0.0015(±0.0010) 




0.90 


0.87 


ax (1+z) 


E/SO 


r\ r\r\'~)C\( I r\ nm o\ 

0.0039(±0.U018) 


-1.6(±0.4) 




0.92 




a+bxz 


S 


0.0069(±0.0025) 


0.0014(±0.0006) 




0.78 


0.70 


ax(l+z) b 


S 


0.0060(±0.0031) 


1.7(±0.5) 




0.92 




a+bxz+cxz 2 


s 


0.0095(±0.0036) 


0.0027(±0.0012) 


-0.0031(±0.0018) 


0.96 








n m = 0.3 




Q A = 0.7 






a+bxz 


E/SO 


0.0049(±0.0026) 


0.0032(±0.0040) 




0.86 


0.83 


ax(l+z) 6 


E/SO 


0.0055(±0.0024) 


-1.1(±0.3) 




0.82 




a+bxz 


S 


0.0071(±0.0022) 


-0.0019(±0.0027) 




0.81 


0.73 


ax(l+z) 6 


s 


0.0043(±0.0025) 


1.5(±0.4) 




0.85 




a+bxz+cxz 2 


s 


0.0063(±0.0042) 


0.0017(±0.0011) 


-0.0012(±0.0015) 


0.87 





p(z) = 0.0039(±0.0018) x (1 + Z y 1 - 6(±0A) . In this case, the 
decrease of elliptical galaxies is ~ 60(±10)%. This behavior 
is in a very good agreement with the prediction from the hi- 
erarchical clustering scenario for this cosmology (Baugh et 
al. 1996; Kauffmann et al. 1996) but differs from the results 
presented in Totani & Yoshii (1998) and Im et al. (2001). 
Interestingly, Daddi (2001) has pointed out that strong dis- 
crepancies in the number density evolution for the EROs 
(Extremely Red Objects^]) can be understood in terms of 
cosmic variance: "it is much probable, on average, to un- 
derestimate the true ERO surface density with small area 
surveys". Maybe a similar explanation also holds for more 
modest redshift E/SO population and this can be of help to 
understand the discrepancies in the number density evolu- 
tion pointed out for different authors. 

For the spirals, the comoving density is p(z) = 
0.0069(±0.0025) + 0.0014(±0.0006) x z, but the KS prob- 
ability is just 0.78. The parameters of the power-law model 
for this family are p(0) = 0.0060(±0.0031) and m=1.7(±0.5) 
with a KS probability of 0.92. Interestingly, for this popu- 
lation a peak in the range z — 0.4 — 0.5 is shown in the 
comoving density obtained from binned data in Fried et al. 
(2001), although they fit only a linear comoving density. 
Probing on this possibility, we have also tested a quadratic 
comoving density: 

p(z) = a + bz + cz 2 (6) 

where the interpretation of these parameters is as follows: 
a=p(0), b=2Ap/z p where z p is the redshift where the co- 
moving density reaches its biggest value and Ap — p(z p ) — 
p(0) and c=-Ap/Zp. Using a quadratic comoving density we 
obtain the highest probability, 0.96, with the next values for 
the parameters: p(z) = 0.0095(±0.0036)+0.0027(±0.0012) x 
z - 0.0031(±0.0018) x z 2 . Notice that this implies a peak of 
the density at z = 0.43. Nevertheless, the value of the co- 
moving density at this peak is just 6% higher that at z = 0. 
Meanwhile the value of the density at z = 0.8 is slightly 
higher (about 1%) than at z = 0. Consequently, contrary 



^ Most of these objects are expected to be E/SOs at redshift 
l<z<2. 



to the E/SOs, brighter spiral galaxies {Mb < —18) seem to 
have a relatively quiet evolution. 

Our values of p(0) for E/SO and spiral galaxies are in 
good agreement with the values that can be obtained by 
using the fit to the Schechter luminosity function (Schechter 
1976) of nearby samples (Marzke et al. 1998). For Mb < 
-18, the local comoving density is p(0)=0.0026 ± 0.0007 
(E/SOs) and p(0)=0.0054 ± 0.0014 (spiral galaxies). 

We have also evaluated the previous quantities assum- 
ing a different cosmology: Q m — 0.3 and J1a = 0.7. In this 
case, our absolute magnitude limit is Mb < —19. We sum- 
marize our results in Table 3. The E/SO comoving density 
at this cosmology seems to no-evolve or slightly decreases. 
This is a similar result to that obtained for this cosmology 
by Totani & Yoshii (1998) and Im et al. (2001) and what 
it is expected from semi-analytical hierarchical models (e.g. 
Kauffmann & Chariot 1998). The results for the spiral galax- 
ies are compatible with no number density evolution. 



6 CONCLUSIONS 

We present quantitative morphology of galaxies in two 
Hawaiian Deep Fields imaged by 6 HST fields. Down to 
the limiting magnitude of our sample, nearly all galaxies 
have spectroscopic redshifts. The morphology has been ob- 
tained by fitting a pure Sersic and Sersic + exponential pro- 
files to the surface brightness distribution of the galaxies. 
Monte-Carlo simulations have been carried out in order to 
determine the limiting magnitude down to which the recov- 
ered structural parameters are reliable. The galaxies have 
been classified according to the B /T ratio. E/SOs systems 
are those with B/T > 0.5. Our simulations suggest an ap- 
parent magnitude limit of 1=23. We have also accurately 
determined the absolute limiting magnitude of our sample 
Mb < —18. The complete subsample is composed by 61 
objects up to z — 0.8. 

The percentage of the different galaxy types in the 
whole sample are in good agreement with those obtained 
in the HDF by visual methods. We have computed the co- 
moving density of the galaxies as function of redshift. For an 
Einstein-de Sitter universe, the comoving density of E/SOs 
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decreases as z increases, in very good agreement with the 
predictions of hierarchical clustering models of galaxy evo- 
lution. The comoving density of spiral galaxies shows a good 
fit to a quadratic form: it grows a ~6% from z — until 
z — 0.43, and then decreases slightly until z — 0.8. This 
fit is compatible with no number evolution. For open or A 
universes, the E/SO galaxies comoving density is compati- 
ble with no number density evolution or a slightly decrease 
as it is expected from semi-analytical models in hierarchi- 
cal clustering scenarios. Density comoving for brighter spiral 
galaxies also remain quite constant at this redshift range. 
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